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Abstract 

We report a theoretical study of the axial Nucleon to Delta(1232) (N A) transition form 
factors up to one-loop order in relativistic baryon chiral perturbation theory. We adopt a formalism 
in which the A couplings obey the spin-3/2 gauge symmetry and, therefore, decouple the unphysical 
spin- 1/2 fields. We compare the results with phenomenological form factors obtained from neutrino 
bubble chamber data and in quark models. 

PACS numbers: 23.40.Bw,12.39.Fe, 14.20.Gk 



1 



I. INTRODUCTION 



The axial A(1232) transition form factors play an important role in neutrino induced 



pion production on the nucleon, in particular at low energies [l, 2, 3, ^, 0]. These form factors 
have been parametrized phenomenologically to fit the ANL 6, 7| and BNL [s, 0] bubble- 
chamber data. In the past, the theoretical descriptions have been done using different 



approaches, for a review, see Ref. 



10|. In recent years, there has been an increasing interest 



on these form factors. They have been calculated, for instance, using the chiral constituent 

121. State of the art calculations within 



quark model and light cone QCD sum rules 

lattice QCD jisl . [l^ have also become available. The possibility to extract the axial — A 
transition form factors using parity-violating electron scattering at Jefferson Lab [15| has 



been studied extensively 
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17|. Present and future neutrino exp eriments could also provide 



further information on these form factors 
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23|. 



Chiral perturbation theory, based on a simultaneous expansion of QCD Green functions 
in powers of the external momenta and of the quark masses, has achieved remarkable success 

nnnn 

in describing the dynamics of the light pseudoscalar mesons at low energies [2J, 1251 . 1261 . |27J . 
The sector with one baryon is more problematic because, as was shown in Ref. |28[], the 
systematic power counting is lost since the nucleon mass is not zero in the chiral limit. 
These problems were first handled in heavy baryon chiral perturbation theory (HB^PT), 



where nucleons are treated semi-relativistically |29|, [SOj. However, in certain cases, this 
approximation leads to convergence problems because the Green functions do not satisfy 
the analytical properties of the fully relativistic theory 3l| . Recently, the systematic power 



counting has also been restored in the relativistic formulation through either the infrared 
or the extended on-mass-shell regularization schemes 



3l| 



33|. 



The explicit inclusion of the A in chiral perturbation theory requires a power counting 
that properly incorporates the A-A^ mass difference, A = — M]\f, which is small com- 
pared to the chiral symmetry breaking scale. Two expansion schemes have been proposed. 
One is the small scale expansion js^ which considers A to be of the same order as the 
other small scales in the theory, i.e., ~ p ~ Z\. The other is the S expansion scheme, 
which counts A differently depending on the energy domain 3^ . Originally, the small scale 
expansion was used in HByPT, while recently it has also been implemented in relativistic 
chiral perturbation theory 36l. l37l|. 
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The vector N ^ A transition form factors, important to understand eA^ (7^) reactions 
and the structure of the nucleon, have been calculated up to next-to-leading order in both 



the small scale expansion HBxPT 38|, |39| and the 6 expansion re 



ativistic baryon xPT 40|, 



4l|. While axial form factors have been addressed in HBxPT 42], no calculation has been 
performe d up to now within the relativistic framework. With lattice QCD results becoming 
available ISj, it is timely to study the axial transition form factors within relativistic chiral 
perturbation theory. 

In this paper, we use the relativistic baryon chiral perturbation theory, including explicitly 
the A resonance, to calculate the axial N ^ A transition form factors up to order 3 in the 5 
expansion. In sect. II, we briefly explain the power counting, the difference between the small 
scale expansion scheme and the 6 expansion scheme, write down the relevant Lagrangians 
up to next-to-next-to- leading order and the appropriate form of the A propagator. Loop 
calculations are performed in sect. III. In sect. IV, we discuss our results in terms of 
the low energy constants and loop functions. In sect. V we compare the results with 
both phenomenological parameterizations and other theoretical calculations. Summary and 
conclusions are given in sect. VI. 

II. POWER COUNTING, EFFECTIVE LAGRANGIANS, AND THE A PROPA- 
GATOR 

A. Power counting 



A fundamental concept of xPT (as Effective Field Theory) is the power counting [24 1. 
It provides a systematic organization of the effective Lagrangians and the corresponding 
loop-diagrams within a perturbative expansion in powers of (p/A^sb)"^^^? where p is a 
small momentum or scale and A^sb, the chiral symmetry breaking scale. In xPT with pions 
and nucleons alone the chiral order of a diagram with L loops, Nt^{Nn) pion (nucleon) 
propagators, and Vk vertices from /cth-order Lagrangians is 

n^PT = 4L - 2N^ _ iV^ + ^ A;Ffc. (1) 

k 

However, in the covariant theory this rule is violated in loops by lower-order analytical 



pieces 



281]. This power counting can be recovered by adopting non-trivial renormalization 
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schemes, where the lower-order power-counting breaking pieces of the loop results are sys 



tematically absorbed into the available counter-terms [31|, 133| . A detailed discussion of the 
renormalization scheme adopted in the present work will be presented together with our 
main results in section IV. 

If the A resonance is explicitly considered, things become more complicated because its 
excitation energy, A = Ma — ~ 0.3 GeV, is small compared to the chiral symmetry 
breaking scale A^sb = 47r/^ ~ 1 GeV. Therefore, there are two small parameters in the 
theory, i.e., 

e = m^/A^SB and 6 = A/A^sb- (2) 
Over the past few years, two different expansion schemes have been proposed, the small 



scale expansion and the 6 expansion. In the small scale expansion 3J], one has m,r 



Z\ ~ p ~ (9(e). In the 5-expansion [35|, to maintain the scale hierarchy <^ A A^sb, 
m-Tr/A^SB is counted as 6^. In this scheme, the power counting depends on the energy domain 
under study: p ~ or p ~ Z\. 

For the study of — * A axial transition form factors in the energy region p ~ Z\, the 
order of a graph with L loops, Vk vertices of dimension k, A^ pion propagators, nucleon 
propagators, A'a Delta propagators, the power-counting index n is given by: 

n = %PT - A^A- (3) 



For a more general discussion, see Ref. 43 1. 



In the present work, we adopt the 6 expansion scheme. As can be seen in the following sec- 
tions, the differences between these two schemes in our case come from vertices proportional 
to m^, which count as 6^ in the 6 expansion and, therefore, have been neglected. 

B. Chiral Lagrangians 

In this section, we write down the relevant A^A^, A^A, and A A Lagrangians and pay 
special attention to the A couplings and the spin-3/2 gauge symmetry. 
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1. Pion-nucleon and pion-pion Lagrangians 



The lowest order pion-nucleon Lagrangian has the following form: 



(4) 



where M^r and qa are the nucleon mass and the axial- vector coupling at the chiral limit, 
is the covariant derivative 

D^N = d^N+[T^,N], (5) 



{u'^idf, - ir^)u + u{d^ - il^)u'] 



and the axial current defined as 



(6) 



(7) 



In the above definitions, = + a^, = — with = t"vJ^/2 and = r°'a^/2 
the external vector and axial currents, where t'^ are the Pauli matrices. The matrix u 
incorporates the pion fields 



u 



Trrll 



U = exp 



TT 



V^TT^ 



V271- -TfO 

with being the pion decay constant in the chiral limit. 

The leading order pion-pion Lagrangian has the following form: 

C^^} = ^Tr [V^UiWUy] + &Ti [xU^ + Ux^ 

where 

V^U = d^U~tr^U + tUl^ 



(8) 
(9) 



(10) 



(11) 



with X = diag(m^, ml). 



2. Nucleon-Delta and Delta-Delta Lagrangians 

The A(1232) is a spin-3/2 resonance and, therefore, its spin content can be described 
in terms of the Rarita-Schwinger (RS) field A^, where /z is the Lorentz index. ^ This field. 



^ We follow Ref. t43i] and write the Lagrangians for the spin-3/2 isospin-3/2 A isobar in terms of the 
Rarita-Schwinger (vector-spinor) isoquartet field = (A++, A+, A°, A~)^, which is connected to the 
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however, contains unphysical spin-1/2 components. They are allowed for the description of 
off-shell Delta's, but the physical results should not depend on them. In order to tackle 
this problem, we follow Refs. [43|, |45|] and adopt the consistent couplings, which are gauge- 



invariant under the transformation 



/\^{x) ^ A^{x)+d^t{x). 



(12) 



A remarkable consequence of the use of the spin-3/2 gauge symmetric couplings is that it 
leads to a natural decoupling of the propagation of the spin-1/2 fields. 

In the following we give the A^A and AA Lagrangians relevant to this work. The lowest 
order Lagrangians in the resonance region are^ 



ihA 



(1) 



2Ma 

Ha 



AA 

where ul = ^Tr (t'^ux) = - j^dxTx" + + ■ ■ ■ , T'^ and T"^ are the isospin 1 /2 to 3/2 and 3/2 
to 3/2 transition matrices, and the totally antisymmetric gamma matrix product as 

given in Appendix A. At second order, there are four terms, i.e.,^ 



2Ml 



(13) 



(14) 



Ma 



ArT'^(9^A^)cj'^'^^ + H.c., 



(15) 



isospurion representation of Ref. 3J| through 



where T" are the isospin 1/2 to 3/2 matrices satisfying T'^T^'^ = 5"^ ~ l/3r°r^ as given in Appendix A. 
With this rule, the on-sheU equivalent form of our consistent couplings can be easily identified with those 



34 



44| 



of Refs. 

^ If one A is put on-shell, the A- A Lagrangian is equivalent to that of Pascalutsa et al. [43|: 



= -^e^'^'"^A,T'^(9,A.)^^ + H.c. 



^ In our study of the axial form factors up to one-loop order the 5^^) and 5^^^ Lagrangians only concern 
on-shell A's. Therefore, they are the same in the consistent coupling scheme of Pascalutsa et al. as those 
conventional Lagrangians in Refs. 3J, |44|. 
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while at third order, there are seven terms^ 



Ma Ma 



3 tTt, 



Ml 



Ml 



(16) 



where uj"^^ = iTr(r'^[L'^., m^]), fl'^" = - d^a""'^. As we will see later, the 5(2) and 5^^) 

low-energy constants (LEG) contribute to the form factors only in particular combinations; 
therefore, the number of independent parameters is smaller than the one appearing in the 
above Lagrangians. 

C. Spin-3/2 propagator 



The most general spin-3/2 free field propagator in D dimensions has the following 



form 



36 



46| 



^ + Ma 
Ml-p^ 
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{D-l) (D-1)(CV-M1^ 



(D-2)(1-C>V 



(17) 



(^-i)(CV-Mi)J 

where C is the spin-3/2 gauge- fixing parameter. In the case of C = 0, the above propagator 
corresponds to the usual Rarita-Schwinger propagator 



i) + MA 

Ml-p^ 



9 



a/3 



{D - 1) {D- 1)Ma 



{D - 2)pV 
{D - l)Ml J 



while in the case of C = oo, it becomes 



^P(rr.\ - + pa/3 . X 



(18) 



(19) 



with the covariant spin-3/2 projection operator defined by 

1 



(^\ _ ^"/3 Tl^ 



KfM = 9^'- 



{D - 1) {D- l)p2 



(i^tV+p"?'^;^) - 



- 4)pV 
{D - l)p2 



(20) 



^ In the small scale expansion scheme, there are two more terms at this order proportional to m^, i.e., 



where x+ and X- ai'C external scalar and pseudoscalar sources. 
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FIG. 1: Feynman diagrams contributing to the axial ^ A transition form factors up to 6^^\ 
The double, solid, and dashed lines correspond to the Delta, nucleon, and pion, while the wiggly 
line denotes the external axial source. 

It should be stressed that due to the spin-3/2 gauge symmetric nature of the consistent 
couplings, our results do not depend on the particular value of the gauge-fixing parameter 



III. THE iV ^ A AXIAL TRANSITION FORM FACTORS 

The N A axial transition form factors can be parameterized in terms of the usually 
called Adler form factors Q : 

+Ci{q')g-^ + ^^q-q^y, (21) 

where A°'^'^ is the third isospin component of the axial current. 

All the diagrams contributing to the ^ A axial transition form factors up to J^^-* are 
displayed in Fig. [H^ Two Kroll-Ruderman like diagrams are not shown since the one with 
an internal nucleon and a AttNN vertex is zero and the other one with an internal A and a 



^ We do not have the diagrams (c), (d), (e), and (f) of Fig. 1 of Ref. 4^, that correspond to tadpole 
diagrams where a pion loop couples to either the ANA (ttNA, Att) vertices, or the pion fields, because 
the contribution of those diagrams are of higher-order in the S expansion scheme. 
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Ak/S./S. vertex contributes as a real constant, which is irrelevant to the present study due to 
the adopted renormahzation scheme. The calculation of the tree-level diagrams [Fig. 1(a)] 
is straightforward: 



(o) 



hA 
- -W9 



Ma 



d. 



1 .J 



Ma 

A 

Ma 



p' ■ qg'^' + hiq^g'^'' - q'^q') + ^^'^"'^ + i^p' ■ q{l ■ qg""" - q'^l') 



Ma' 



1 ■ q{p' ■ qg"'^ - q^'p'^) + n— 7 • qp' ■ qg"'^ + tttp' ' q^' ' qs"'^ ~ q"p''^) 



Ma 



Mr 



+j^p'-qp'-qg"^ 



(22) 



where p', p, and q are the momenta of the A, the nucleon, and the external source. We 
assume that both the external nucleon and A are on-shell, which yields p' • q/M^ — (M^ — 
M|r + g^)/(2MA) ^ A and 7 • g = Ma - Mjv = A, where we have neglected the and 
terms which, strictly speaking, are of higher order than the chiral order of the corresponding 
Lagrangian. 

In the following wc explicitly show how to calculate the loop diagrams: 
Diagram Fig. 1(c) reads 



hAg: 



i 1 



(87r/,)2MA 



with 



iG^^. = (27r/.) 



4-D 



K7'"'^c [f-H Mn] Yi5 l^-H Mn] H 



(c) J -^2 [^2 -rril + ie] [{p - kf - + ie] [{p' - kf - + ie] 

where //, the renormahzation scale, is set to be Ma. 



(23) 



(24) 



Diagram Fig. 1(e) reads 



gAhAHA 1 
, (87rA)2 Ml 



with 



zG;'; = (27r/.) 



4-D 



d^k ie'^''^%kcSad{p' - k)Y'^''{p' -k)e{i>-lj^ + Mjv)^75 



- m2 + ie] [{p - ky - Mj, + ie] 



(25) 



(26) 



Diagram Fig. 1(g) reads 



4«M,3 _ 1, /2 
^(9) -6V3 



inU)^Ml 



(27) 



with 

' ) ^ J Z7r2 [A;2 - m2 + ze] [{p' - kf - Mj, + ze] 



(28) 



Diagram Fig. l(i) reads 



^« "9V3 



hAH\ 



(29) 



with 



dPk 



X 



, ■ . (30) 

jp' - k)pk^Sa,{p' - k)^^'^Yl'llcn{p' - k)m{p - kYS'\p - fc)7/dg(p ' k)f k^ 

— ml- + iel 



[A;2 — ml + ie] 

In the above equations, S^^{p) is the spin-3/2 propagator defined in Eq. f|T71) . Since the 
couphngs we used are spin-3/2 gauge symmetric, our resuhs do not depend on the specific 
value of the gauge fixing parameter. 

These loop functions are quite complicated, particularly the ones including A internal 
lines. In practice, we adopt the conventional Feynman parametrization method (see Ap- 
pendix B) and calculate these loop functions numerically. The manipulation of the Dirac 
algebra has been performed independently with FORM 48|| and FeynCalc [49]. The resulting 
Feynman parameter integrals are listed in Appendix C. Whenever possible, the numerical 
results have been checked using the FF library 50| through the LoopTools interface 51|. 

The one-loop results contain only four different Lorentz structures (due to the constraints 
Aq,7" = and Aqp'" = 0), i.e., 7^^^", g""^ , and In accordance with the Adler 

formulation of Eq. (1211) . we can identify the corresponding Lorentz structures and group the 
results as 



(c) "I" (e) (g) (i) 



+g5{q')9'"' + 96{q')q"q'' 



qyn 



(31) 



couplings: m^. Mat, Ma, /tt, gA, hA, and Ha- Here, we adopt the following values: 



It is interesting to note that these loop results depend only on known masses and 
Ma 

= 0.139 GeV, Mn = 0.939 GeV, Ma = 1.232 GeV, = 0.0924 GeV, gA = 1-267, 
Ha = 2.85, and Ha = {9/5)gA- The value of Ha is obtained from large Nc relations and its 
uncertainty is discussed below. In other words, the dependence of the loop functions are 
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genuine predictions of the present work, in contrast with the S^^^ and S^^^ tree level diagrams, 
which contain basically unknown low energy constants: di, d2, d^, d^, fi, /2, /s, f^, /s, /e, 
and fj. Some of these LEG, (c/3, ^4, /s, /?), also appear in pion-nucleon scattering and 



44|. 



could, in principle, be extracted from there 

Apart from diagrams (a), (c), (e), (g), and (i), the external axial source can also couple 
to a pion and interact through it with the system. These are the so-called pion pole terms 
(diagrams (b), (d), (f), (h), and (j)) and are calculated below. 

The Lagrangian responsible for the coupling of the external axial source with the pion at 
second order is 

>C(') = -Ud.n'^a'''''. (32) 

With this and the low-energy counter terms given above, we can easily write down the 
pion-pole contributions: 



pion— pole 



2 



3g2 



mi 



h 



A . di , 2 /s / 



■^{p' -qY - {95 + gaq'^) 



(33) 



with and Qq the loop functions calculated above. 



IV. RESULTS AND DISCUSSIONS 



In this section, we present our results for the form factors in terms of the LEG and 
the loop functions g^, g^, g^, and gQ (Table [T]). It should be mentioned that the Partially 
Gonserved Vector Gurrent (PGAG) relation 

Ci + ^q'U^o = (34) 

holds up to every order in our xPT study, which can be easily checked from Table [B 

As mentioned above, the one-loop results are free of unknown couplings, but the LEG 
are basically not known. Since these LEG always appear in particular combinations, we can 
introduce di = di — {f^ + /e)^, ^2 = ^2 — /a^, and ^3 = ^3 -|- c?4 — (/s -|- /7)Z\ and treat 
them as free parameters. Therefore, effectively, we have five unknown constants: di, ^2, d^, 
/i, and /2. 

From Table [H we can conclude that 
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TABLE I: The axial transition form factors in relativistic baryon chiral perturbation theory; di, 
d2, ds, d^ are order 2 LEG (in units of GeV^^) while fi, f2, /s, fi, /s, fe, h are order 3 LEG 
(in units of GeV~^); gsiq'^), g^iq'^), 95(9^), and g&{q^) are the one-loop contributions as defined by 

Eq. dan). 



FF 



5(1) 



<5(2) 



5(3) 



V 2 Mjv 





-d2 


V2 M-j, 





-di/MA 




Ha 
2 


-(d3 + d4)Z\ 


\J2 Ml 


/ia/2 





/3Z\ + 53((?') 

(/4 + /6)^/MA+54(g') 
(/5 + /7)Z^2 + (^^+^^y+^^(^2) 



T" 



(a) At order 5(1), Cg^ = 0, Cf = 0, and 



2hA ~ 
3 2 ~ 



1.16 with hA = 2.85 from Ref. |46|, 



which is determined from the A-resonance width, Fa = 0.115 GeV. Furthermore, Cq 
is related to through the pion-pole mechanism, i.e.. 



= a 



5 9 9 * 

— 



(35) 



(b) At order 5^'^\ C^, C^, and receive a finite constant contribution. The above 
relation, Eq. ( 135|) . between and Cq still holds. 

(c) At order 5(3), the LEG give constant contributions to all form factors, and dependent 
contributions to and Cq. The one-loop diagrams start at this order. 

Before presenting the loop results we specify our regularization procedure due to the 
complications with the power counting mentioned in Section II. A. The loops are regularized 



in the MS scheme, subtracting in addition the real part of the contribution to the form 
factors at q^=0. Since there is no counter terms linear in at 5^'^\ this procedure guarantees 
to recover the power counting in all form factors. 

We show in Fig. [2] the one-loop contributions to the form factors C^, C^, C^, and 
Cq (except the pion-pole diagrams which only contribute to Cq). One can see that only 
diagrams c, d (N-N) and g, h (A-A^) from Fig. 1 contribute to the imaginary part of the form 
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FIG. 2: (Color online) One-loop contributions to the form factors C^, C^, C^, and Cq. The pion- 
pole diagrams, which only contribute to Cq, have not been included. The N-N, N-A, A-N, A-A 
labels denote the contributions of diagrams with nucleon-nucleon, nucleon-Delta, Delta-nucleon, 
Delta-Delta internal lines. 

factors, with A^-A^ being dominant. One also finds that and receive relatively small 
corrections from the one-loop calculation, whereas gets a relatively large one coming from 
the A-A diagrams (diagrams i, j). This seemingly large dependence, however, suffers 
from the uncertainty related to the vrAA coupling Ha because the A-A loop contribution 
is proportional to H\. 

In Fig. [31 the loop contributions from all diagrams to each form factor are added. Clearly, 



13 





-q^ KGeV/cf] 



-q^ [(GeV/c)^] 



FIG. 3: (Color online) One-loop contributions to the form factors C^, C^, C^, and Cq. The 
pion-pole diagrams, which only contribute to Cq, have not been included 

one can see that has the largest imaginary part; the second; next is the C^, and Cq 
receives the smallest contribution. 

Without the one-loop contributions, Cq can be easily separated into a non-pole part and 
a pion-pole part, i.e.. 



1 



q2 _ ^2 g 



with 



2 



6 



d 



■^(/i + /2)^6-^log(C5^)|,.=o. 
9-kNA , dq 



This is equivalent to the HBxPT result of Ref. [42l] 



6 1 



NA 



h + h K 

2 M, 



+ C2 



with the correspondence Qt^na = Q-kNA and (/i + 



N 



2 Afjv 



+ C2 



(36) 

(37) 
(38) 

(39) 



V. COMPARISON WITH OTHER APPROACHES 
A. Phenomenological fits 

Bubble chamber neutrino data have been used to extract information about the axial 
N ^ A form factors {4,9, 



52 



53| . However, there are some important limitations. First, the 



cross section is basically dominated by the form factor and shows very little sensitivity to 



14 



C;f4 6- Second, the statistics is quite low and, furthermore, the two available data sets from 
BNL jsj] and ANL [3] are clearly different. Finally, it is difficult to disentangle the A from 
other background pion production processes 3|, |J]. Therefore, all these works make some 
additional assumptions. A set of them often found in the literature^ is: C^=0, = —jC^, 
and Cq is related to through Eq. (l35ll . In this way only is fitted to the experiment. 



.5 tiiiuug,!! 1:74. vbJi^- J-ii tiiiD vva_y uiii_y u/5 

As an example, we can take Kitagaki et al. [9'| where it has the following functional form: 



2 



q2 



2 



Ml 



-2 



(40) 



with ^^(0) = 1.2, as = -1.21, h = 2 GeV^, and Ma is fitted to data yielding Ma = 
1.28]^Q'j'o GeV. We will refer to this set of form factors as Kitagaki- Adler (KA) form factors. 

As we have shown above, there are 5 independent parameters in the 6 expansion scheme 
up to chiral order 3: di, d2, d^, fi and f2- We fix them in such a way that the real part of 
our form factors reproduces many of the features of the KA ones. To obtain = 0, we 
set (^2 = 0; therefore its contribution comes only from loop calculations which are of chiral 
order 3. Strictly speaking , Eq. fl35|) is not fulfilled at order 5^^^ but, if one neglects the 
small loop contributions, it can be satisfied by taking /i = 0. Correspondingly, the relation 
C^(0) = -lC^{0) fixes di; C^{0) = 1.2 fixes 4- The only LEG left, /a, is then fixed to 
reproduce at = 0. 

The results obtained this way are shown in Fig. HJ with the following parameter values 
di = -0.514 GeV-\ rfs = 0, 4 = 0.153 GeV-^, /i = 0, and = -2.184 GeV'^. One 
can clearly see that the calculated and Cq are in good agreement with the KA form 
factors. On the other hand, the dependence of is much weaker that the one assumed 
m KA, = and we cannot accommodate their results at order 6'^^\ For C^, the 

dependence is also very weak (compared to C^) . In Fig. HI the dark shadowed area indicates 
a modification of Ma within its uncertainties as given in Ref. qI]. As we mentioned above, 
the dependence on is rather sensitive to the coupling constant Ha- This can be easily 
seen from the light shadowed area in the upper-left panel of Fig. HI which covers the region 
of Qa < Ha < {9/5)gA- The form factors, C^sg, on the other hand, are less sensitive to the 
value of Ha- 



° This choice originates from the analysis of Refs. of Adler's results obtained using dispersion 

relations [551 . 
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FIG. 4: (Color online) Comparison with the Kitagaki-Adler form factors. The dark shadowed area 
indicates the uncertainty of Ma = 1. 28^^050 GeV determined in Ref. jsj]. The light shadowed 
area indicates the sensitivity of the results to the vrAA coupling Ha, which covers Ha = (O/S)^^ 
to Ha = QA- 

A word of caution is in place about the comparison of and with the KA form 
factors. In xPT, the leading order counter terms linear in contributing to and 
appear at chiral order 4. A fair comparison with the phenomenological fits (particularly 
the dependence) should, in principle, be done at order 4. However, the 5*^^^ xPT results 
might give us a clue on the magnitude of the q^ dependence of and Cf. Furthermore, 
if we believe in the phenomenological assumption, or the results of other approaches, the 
difference between the third order xPT results and the results of other approaches might 
help us estimate the value of the corresponding fourth order LEG. Indeed, the upper panels 
of Fig. m indicate that small 5^^'> corrections hiq^ with natural values for the LEG hi can 
reproduce the slope assumed for and C/ by the KA ansatz. 

we a.o „ot,ce that a .cent analys. U obta,ned a s,„a„e. value fo. bv i„e,„di„g 

non A contributions and fitting to the low invariant mass ANL data. In the present xPT 
study, we do have the higher-order contributions, d^, which could alter C^(0) within such 



16 



a range; however, the same LEG appear in pion-nucleon scattering processes. A combined 
analysis is mandatory to determine whether one can accommodate the small C^(0) obtained 
in, for instance, Ref. |^. This is left for future studies. 



B. Quark models 

There have been many studies of the N A axial transition form factors in various quark 
models, both relativistic and non-relativistic. For a brief review of quark model studies, we 



refer the readers to Refs. 



3, 



Compared to dynamical model studies, a feature of most 
quark model calculations is that the obtained form factors are real due to time-reversal 
symmetry, while in dynamical models, like our xPT study, these form factors are in general 
complex due to the opening of the pion-nucleon channel. 

Quark model results are in fact quite scattered. Taking, for instance, the models discussed 
in Ref. [HI], we observed that the prediction of C^(0) runs from 0.81 to 1.53, C^(0) runs 



from —0.66 to 0.14 and C^{0) runs from to 0.05. These models also obtain the non-pole 
part of Cq whose value at = ranges from —0.72 to 1.13. We could use these results 
to extract our constants although the large differences between them do not allow to reach 
solid conclusions about their values. From C^(0) one gets d^, and from its slope dC^/dq"^ at 
= 0, (/i + /2). This fixes the non-pole part of C^(0) (neglecting the one-loop corrections) 
since 

C^(-^--^^^^\0) ^ ^MUfi + f2) , (41) 
which is nothing but a direct consequence of PCAC. Using the quark model calculation 



of Ref. 



11] for we obtain (J^(^°^-p°^'^) (^q^ ~ _2. This value is almost a factor three 
larger in magnitude than the one obtained directly from that model in spite of the fact that 
it implements PCAC at the quark level by introducing one- and two-body axial exchange 
currents. 

Analogously, we can use quark model results for and at = to obtain di and 
d2. The smallness of C^{0) predicted by all calculations points towards a d2 close to zero, 
in agreement with the phenomenological assumption. The situation is much more uncertain 
with di, both in sign and magnitude. In Fig. (H the dependence of the real parts of 
and in our calculation, which at order 5^^^ is dictated by the loops, is compared to several 
quark models. As in the case of the KA form factors discussed above, we can expect from 
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FIG. 5: (Color online) Comparison with the non relativistic Isgur-Karl and D-mixing quark model 
results of Ref. and those of the chiral quark models of Refs. [Hi] and 56 1. 



this comparison that next order terms linear in with small {natural) values of the LEG 
are sufficient to eliminate the discrepancies in the low behavior with any of these quark 
models. 



C. Lattice QCD results 



Recently, the N —>■ A axial transition form factors have been studied in lattice QCD 



Some major conclusions are (i) and are suppressed compared to and Cq, 
and (ii) can be described by a dipole ansatz C^(0)/(1 + Q^/M^)^ but with a smaller 
C^(0) and a larger Ma (> 1.5 GeV), compared to the Kitagaki-Adler form factors. These 
results should be taken with caution because of the still relatively large pion mass (> 350 
MeV) used in the study. 

In principle, xPT is the perfect tool to extrapolate the lattice QCD results to the physical 
region. Meanwhile, one can also fix the unknown couplings to the lattice QCD results. Due 
to the regularization method we used and the fact that the lattice data points are still scarce, 
we will leave this subject to the future. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the axial — > A transition form factors up to one-loop order in relativis- 
tic baryon chiral perturbation theory with the 6 expansion scheme. The adopted Lagrangians 
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including the A(1232) are consistent, i.e., spin-3/2 gauge symmetric, which automatically 
decouples unphysical spin-1/2 fields. Consequently, our results do not depend on the specific 
value of the gauge-fixing parameter that is present in the most general spin-3/2 propagator, 
and avoid various problems related to inconsistent couplings. 

The form factor exhibits the richest structure in our study. It receives contributions 
starting at chiral order 1, at which we find that C^(0) = -^1^ ~ 1-16 for Ha — 2.85. At 
higher orders, this value is modified by low energy constants that are unknown but which 
also appear in pion-nucleon scattering. At chiral order 3, this form factor gets (f dependent 
contributions, some of them complex. Actually, we find that has the largest imaginary 
part among the four form factors. We also obtain that, up to chiral order 2, Cg = ^2 !lq2 ■ 
At order 3, Cq has a non-pole contribution whose value at = is related to the slope of 
at = 0. Assuming natural values for the LEG, this non-pole part is small compared 
to the dominant pion-pole mechanism. 

Both and Cf start at chiral order 2 and get their dependence at order 3 from the 
loops. For C^, we find a small dependence, which is quite sensitive to the ttAA couphng 
constant. Ha- On the other hand, its imaginary part, coming mainly from the A'^-A'^ internal 
diagram, is finite (~ 0.03 at = 0) and has a mild dependence. This suggests that is 
small (compared to (^4^5 g) but not necessarily zero. The dependence on q^ is also found 
to be rather mild at order S'^^\ 

We have compared our results with a phenomenological set of form factors used in the 
analysis of neutrino-induced pion production data and also with different quark model cal- 
culations. They could be used to extract the low energy constants but the scarcity of data 
and the large differences between quark model results make it difficult to come to solid con- 
clusions. In the case of and specially Cf, the comparison should, in principle, be done 
at order 4 where corresponding leading order counter terms linear in appear. Neverthe- 
less we can say that reasonable agreement with all these approaches can be obtained with 
natural values of the LEG. 

Future experiments with electron and neutrino beams, combined with the analysis of 
pion-nucleon scattering data, can shed more light on these form factors. The extrapolation 
of lattice QGD results to the physical region should also be pursued. 
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Note added in proof: After submitting this paper, a new preprint [57] appeared that 
studies the N ^ A axial form factors up to one-loop order in HBGhPT using the small scale 
expansion scheme. Within this framework, there is no dependence coming from the loop- 
functions. This supports the smooth dependences found in the present work. Namely, 
the dependence of the loops in our relativistic framework is counted as of higher-order in 
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VIII. APPENDIX 

A. Isospin transition matrices and antisymmetric Gamma matrix products 

The isospin 1/2 to 3/2 and 3/2 to 3/2 transition matrices and T" appearing in the 
A'^A and AA Lagrangians are given by: 
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The totally antisymmetric Gamma matrix products appearing in the consistent NA and 
AA Lagrangians are defined as: 

(48) 



7^^ = 2 [7^71, 



with the following conventions: g'^'^ 
^7o7i7273- 



-Z£'^^'"^7<x75, 



diag(l, -1,-1,-1), £0,1,2,3 = 



(49) 
(50) 

1, 75 = 



B. Loop functions 



In the calculation of the loop diagrams, we have used the following (i-dimensional integrals 
in Minkowski space: 



/ 



d/2 



T{X-n + e-2) (-l)"^?^- 



•«2n 



(51) 



■(A^2_p)A 2"r(A) (A<2)A-n+.-2 

with (^^i-- "2n — gaia2 _ _ gOi2n~iOi2n _|_ _ combiuatiou symmetrical with respect to the per- 
mutation of any pair of indices (with (2n — 1)!! terms in the sum 
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The M"^ that appear in the calculation of the N-N, N-A, A-N, and A-A internal dia- 
grams are, respectively, 



^xml-x{l-x- y)Ml + {1 - x - xy)M% - y{l - x - y)q^ - ie, (52) 

^xml + {l-x){l-x- y)Ml + y{l - x)Ml - y{l - x - y)q^ - ie, (53) 

•^AAT — xml + (x'^ + xy — X + y)M\ + {I — x — y — xy)M'^ — y{l — x — y)q^ — ie, (54) 

^xml + {l-2x + x^ + xy)Ml - xyM% - y{l - x - y)q^ - ie, (55) 



where x and y are Feynman parameters. 

C. Feynman parameterization integrals 

We present below the loop integrals, diagrams (c), (e), (g) and (i) of Fig. 1, cast in the 
Feynman parameterization. We use the following notation: Q ^ is the M5'-regularized 
contribution of the loop to Cf with X and Y being the baryons in the internal line (in 



this order), M\y=^'xy / r = M^/Ma, /i^ = m^/MA, and Q^"" = Q^'^/M^ (with 



The couplings are contained in the constants Cxy- 
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a 



Then, the expressions of the loop functions are: 

{NN) _ CnN 
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+ l)a; + y (y2 + y - 2))) j -1^ - p + ((l7 - GQ^) y - 3) a;^ + (y {2Q\Q - lly) 
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- 1) + y - - 5) + (y - 4)y + 1) - 1) r + 3x^ + x^ ((5^(2 - Ay) + 5y - l) 
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+ 3QS{x + y- l){y{lAy - 13) + x{lAy - 9) + 3) + Q' (2(7 - 3Ay)x^ + 2 {y (3(14y - 9y 

- 2{y + l)r - 82y + 35) + 2)x'+ {y {6{y{19y - 26) + 9y + {y{lAy + 23) - 20)r + 4(11 

- 31y)y + 53) - 26)x + y {3(y - l)(2y - l)(5y - 3)r' + {y{y{18y + 11) - 40) + 17)r 
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+ 7/(57-47/(7^ + 3)) -29) + 8) -4) 
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(-16a; - 62y + 12Q'{y{13y - 12) 
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X ^\^(r' + 2{Q'-l)r'+{Q' + l)')y{x + y-l){x'+{r' + Q'-2{r' + Q'-l)y + l)x' 

+ {(y - l)yr^ + {-2y^ + 2Q\y - l)y + l) - 2yr + r + y (Q^ - 1) - 4y + 2) + y + 3) 

- 1 ) + {y{y{2y - 3) + 1)Q^ + {y {{y{2y - 3) + ly + 2{y - l)r - y{2y + 1) + 2) - l) 
+ (r - l)((r + l){2r{y - 1) - 2y + l)y + 1)) x + Q\y - \)y {{Q\y - 1) - 2) y + r{2y 

- 1))) 1 \ (lla^V^ + 8xV^ - 9xr^ - 6xV + 20xr^ - Ar^ + lOxV^ + QxV^ - 12x^2 

- xr^ + 4r^ - lOxV - 4xV + 32a;V - 12xr + 4r + lOx'' + lOx^ + 2(r - l f{r + l)^((15r 

+ 16)x + 15)y^ - ISx^ + (r - lf{r + 1) (6x(7a; - 8)r^ + (30a;2 + x - 42) + 3 (-20x2 + x 

- 4) r - 54^2 + 4x + 31) + 2x - (r - 1) (9x(3x - 2)r^ + 2(x(10x + 9) - 6)/ + {x{x{52x 
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+ r{-12y + 9r{Uy - 9) + 7) - 7) + 2)x2 + (y (6(lly(3y - 4) + 15)r2 + ((71 - 24y)y - 38)r 
+ {y- l){96y - 113)) - 17) x + y {l8{y - l)(r - 2ryf + y{y{7Ay - 191) + 147) + 2r(ll 

- 2y(y(3y - 16) + 17)) - 34) + 4) + g2 (3^ (3(14^ _ g)^;^ + 4(^(12^ _ 17) + 6)x + (y - l)(2y 

- l)(5y - 3)) + 2y ((7 - 12y)x^ + ((60 - \\y)y - 37)x + y{y{y + 39) - 54) + 17) + (2(11 

- h2y)x^ + 2((22 - %\y)y + 5)x2 - {y{y{\2^y + 23) - 76) + 26)x + |/(i/((33 - ^2y)y + 21) 

- 14) + 4) r2 + ((22y - 6)x^ + (42^^ - 347/ + 2) x^ + {y{2y{9y + 32) - 79) + 26)x - (y - 1) x 
X {y{2{y - 45)?/ + 51) - 12)) r + (x + y - 1) (32y^ - 4(13x + 27)y^ - x(74x + 41)y + 22y 
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+ r(-28x + r(32x + 6r(13x - 6) + 3) + 6) + 67)y - 12Q^{y{13y - 12) + x{13y - 4) + 3) 
+ g2 (_6 (g(i _ 27/)2 + 4x(13y - 4)) r2 + 2 (46y2 - 66y + x{A6y - 13) + 22) r - 226y^ 
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